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Abstract. The influence of controlable parameters like temperature and wavelegth on the trajectories 
of light in a nematic liquid crystal with topological defects is studied through a geometric model. The 
model incorporates phenomenological details as how the refractive indices depend on such parameters. 
The deflection of light by the topological defect is then shown to be greater at lower temperatures and 
shorter wavelengths. 

PACS. 61.30.Jf Defects in liquid crystals - 61.30. Dk Continuum models and theories of liquid crystal 
structure - 42.70.Df Liquid crystals - 78.20.Fm Birefringence 



1 Introduction 

Topological defects appear everywhere in physics. Rang- 
ing from gravitation and cosmology to Bose-Einstein con- 
densates they are associated to symmetry-breaking phase 
transitions. In particular, they appear in the isotropic- 
to-nematic phase transition of calamitic liquid crystals 
in the form of hedgehogs, disclinations, domain walls or 
more complicated textures [TJ[2]. The deflection of light 
rays by these defects is an old issue. Grandjean [3J, al- 
ready in 1919, calculated the light paths of extraordinary 
rays passing by disclination lines [4] . More recently, Joets 
and Ribotta [5] used a geometrical model to describe the 
light propagation in anisotropic and inhomogeneous me- 
dia like liquid crystals. Inspired by their work we studied 
the propagation of light near disclination lines in nemat- 
ics from a geometric point of view @I2IE]- In particular, 
in we observed lensing effects due the deflection of the 
beams in regions near the defect. There, we calculated the 
light trajectories from a geometry resulting from the ap- 
plication of Fermat's principlcQ associated to an effective 
refractive index N [10] given in terms of n and n e , the or- 
dinary and extraordinary refractive index, respectively. In 
a few words, what we did was to consider the bent light 
rays as geodesies of a model space of unknown geome- 
try. By identifying Fermat's principle with the geodesic 
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1 Kline and Kay [9] were probably the first to prove that 
light rays in inhomogeneous anisotropic media are extremals 
of Fermat's functional. 



variational principle we were able to find the effective ge- 
ometries for each defect studied. Knowing the effective 
geometry, the geodesies were obtained numerically. Inci- 
dentally, a geometric model describing elastic properties 
of nematic liquid crystals which also leads to an effective 
geometry, appeared recently in the literature [13j . 

Since the refractive indices of a nematic depend both 
on the temperature and on the wavelenght of the light, a 
more realistic model incorporating these effects is in order. 
Li, Gauza and Wu [H] modeled the temperature effect on 
the nematic refractive index based on Vuks equation |15j 
and, by fitting their final expression to experimental data 
of selected materials, found the unknown coefficients. The 
same group derived Cauchy formulae (see page 100 of [TO] ) 
for the refractive indices of a nematic sample as a function 
of the wavelenth of the light [16] . Again, by fitting their 
final expression to experimental data, all coeficients were 
determined. In this work, we incorporate their models to 
our geometric model for propagation of light in nemat- 
ics with topological defects, in order to study tempera- 
ture and wavelength effects. Although in [6] we studied 
both the symmetric (k = 1) and asymmetric (k 1) de- 
fects, without loss of generality, we keep our analysis here 
mostly for the symmetric cases since their effective geom- 
etry is simpler and more intuitive than the asymmetric 
cases. Since, in all cases, the effective geometry depends 
only on the ratio a — n e /n a it is how the temperature and 
the wavelength affect this ratio what matters. In section 
II we review the geometric model and discuss the simplest 
effective geometry for light traveling by a disclination. In 
section III we use the results of P3] and [TB] to show how 
a is affected by temperature and wavelength and study 
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the effect of the variation of this ratio on the light paths 
near selected defects. 



2 Geometric Model 

Disclinations in nematics are classified according to the 
topological index (or strength) k which gives a measure 
of how much the director rotates as one goes around the 
defect. That is, the director configurations, in the plane 
x — y, are given by pp 



<p(0) = k9 



(1) 



where (p is the angle between the molecular axis and the 
x-axis, 9 is the angular polar coordinate and c = <p(0). 
Selected director configurations can be seen on Figure 11.4 
of pp. We assume the disclinations are straight and lie 
along the z-axis and the light rays propagate in the x-y 
plane so, effectively, we have a two-dimensional problem. 

We consider an optical medium constituted by a ne- 
matic liquid crystal with disclinations [2] , where the effec- 
tive geometry for the light is defined by the line element 
(equation (25) of [5]) 



ds 2 — [n cos [(k — 1)4> + c] + n e sin [(k — l)cf> + c]} dr 
+ {n 2 sin 2 [(fc — 1)4> + c] + n 2 cos 2 [(fc — l)<j> + c]} r 2 d(f> 



{2( 



2 



n e — n o) sin[(fc — l)<f> + c] cos[(fc — 1)4> + c]} rdrdcj>. 

(2) 

The metric ([2]) was obtained by identifying Format's prin- 
ciple with the variational principle that determines the 
geodesies in Riemannian geometry. Let 



T 



Nd£, 



(3) 



where, dt is the element of arc length along the path be- 
tween points A and B and the effective refractive index 



N 



n 2 a cos 2 (3 + n 2 sin 2 /?, 



(4) 



where (3 = (n, S) is the local angle between the director n 
and the Poynting vector S. Then, among all possible paths 
between the generic points A and B, Fermat's principle for 
the extraordinary rays grants us that the path actually 
followed by the energy is the one that minimizes T . 

In Riemannian geometry the line element ds depends 
on the position coordinates x l of the point of the manifold 
under consideration. That is, 



ds 2 = gjjdx l dx\ 



(5) 



where gij = gij{x l ) is the metric tensor. The geodesic 
joinning points A and B in such manifold is obtained by 
minimizing ds, just like Fermat's principle. This leads 
to a nice interpretation of the light paths as geodesies in 
an effective geometry [10] . Thus, we may identify 



N 2 d£ 2 = 



,9a 



dx l dx 3 . 



(0) 



The meaning of this equation is the following: the line 
element of the optical path, in an Euclidean space with 
refractive properties, is identified with the line element of 
an effective geometry characterized by . 

In [7] we showed that the effective geometry for the 
vortex-like k = 1, c = 5 disclination is that of a cone. The 
effective metric for this case is obtained by substituting 
these values in metric @ and rescaling the coordinate r 
to p = n e r. That is, the two-dimensional line element for 
this effective geometry, in polar coordinates, is [6] 



ds 2 = dp 2 



a 2 p 2 d6 2 , 



(7) 



where a = n a /n e is the ratio between the refractive in- 
dices. The geodesic equation in a Riemannian space like 
the cone is [IT] 



d 2 x l 
dt 2 



\ - ni dxi dx k 



dt 



0, 



(8) 



where t is a parameter along the geodesic and r l - k are the 
Christoffel symbols, given by 



dgkm , dg 
dxJ 



'"J 



9 9. 



3k 



dx k dx m 



(9) 



For metric (0 equation ([8} reduces to the coupled system 
of ordinary differential eaquations 



d 2 p 2 (dO 
^- aP {di 



(10) 



and 



dH_ 

dP 



0. 



(11) 



2dpd9 
p dt dt 

The solution to the coupled system (fTO]) and (fTT|l is 
easily obtained [IT] : 



Pit) 



C 2 
Ea 2 



2E{t + DY 



n/ . 1 f2Ea(t + D) 

9{t) = - arctan i '- 

a \ c 



F 

a 



(12) 



(13) 



where C, D, E and F are integration constants. 

In figure 1 we show the light paths in the nematic 
medium with the k = 1, c = ? disclination as given by 
(Ti"2"]) and i(T5|) . In figure 2 the geodesies on a cone are shown 
for comparison. 

Metric ([7J describes a cone. Fig. 3 shows the making 
of a cone from a planar sheet where an angular section 
was removed with posterior identification of the edges. 
If 7 is the angle that defines the removed section then 
the remaining surface corresponds to an angular sector of 
2na = 2tt — 7. This is exactly what metric ([Jj describes. 
The incorporation of the term a 2 to the Euclidean metric 
in polar coordinates makes the total angle on the surface 
be Jg 27r ad9 = 2na < 2tt, since n Q < n e . It is clear then 
that a tells how "pointed" is the cone. The closer a gets 
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Fig. 1. Light trajectories in a nematic liquid crystal with a 
topological defect given by a disclination k = 1 and c = tt/2. 



3 Refractive index variation 

The refractive indices n and n e of a nematic liquid crystal 
depend both on the temperature (T) and on the wavelenth 
(A) of the light. In this section, based in [14] and [16], we 
analyse how these parameters affect the ratio a ~ n /n e , 
which characterizes the effective geometry associated to 
disclinations. By changing either T or A, a is changed and 
so is the effective geometry. This causes a deformation of 
the geodesies associated to the light rays in our model. 

In [14j we can find expressions to the ordinary and 
extraordinary refractive index given in terms of the bire- 
fringence An and of its average value (n), such that 



n a = (n) - i/in, 



(17) 




Fig. 2. Geodesies on the cone. 




Fig. 3. Conical surface of angular deficit 7. 



to 1 the flatter is the cone. For a = 1 the cone turns into 
a plane. 

The solutions ([12]) and ([13]) still hold for the radial 
defect with k = 1 and c = since, in this case, the line 
element @ reduces to 



ds 2 = dp 2 + -\p 2 d0 2 , 
or 



(14) 



where, a = n a /n e still. Consequently equations (|12[) and 
([T^]) become 



P(t) 



C 2 a 2 
E 



2E(t + D) 2 , 



f2E(t + D)\ r . 

6(t) = a arctan ) + Fa, 

V ca 



(15) 



(16) 



(n) + -An. 



(18) 



In ([T7]) and ([TS]) the behavior of (n) as function of the tem- 
perature [14] is given through a linear dependence given 

by 

(n) = A - BT, (19) 

where the parameters A and B are obtained experimen- 
tally. 

The birefringence can be written in terms of the ap- 
proximated [18j order parameter S = ( 1 — ijr J as 

An = (An) U-^ r y , (20) 

where (An)o is the birefringence at T = OK, (3 is a con- 
stant associated to the material and T c is the isotropic- 
nematic transition temperature. 

Therefore, substituting the equations ([19]) and ([20]) 
into (|T7[) and ([T5]) . we have 



= A-BT 



A- BT 



(An) 



2{An) 



1 - 







T 



(21) 



(22) 



The liquid crystal considered was the 5CB or 4-cyano- 
4-n-pentylbiphenyl and the wavelength of the incident beam 
was 589nm [M] . For this material the parameters are given 
in the table below obtained from [14]. The parameters A, 
(3 and (An)o are adimensional. 



A 


B 





(An) a 


T c 


1.7546 


0.0005360 K- 1 


0.2391 


0.3768 


306.6 K 



where, as before, C, D, E and F are integration constants. 



In figure 4 we show the ratio a = n a /n e as a function 
of the temperature as obtained from equations ([2"T]) and 
(|2"2"]) . As shown in [7] these parameters are associated to 
the effective geometry perceived by the light traveling in 
the vicinity of k = 1 disclinations. For these defects the 
geometry is conical with the radial disclination (fc = 1, 
c = 0) behaving as a negative curvature cone and the 
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Fig. 4. Effective geometry parameter a as function of the tem- 
perature for 5CB in the nematic phase at 589 nm. 



Fig. 5. Effective geometry parameter a as function of the 
wavelength for 5CB in the nematic phase at 25.1 °C. 



vortex- like disclination (fc = 1, c = 7r/2) like the ordinary 
cone. The value a = 1, reached at T c corresponds to the 
Euclidean geometry which describes the isotropic phase. 

Next, we consider the wavelength dependence of the 
effective geometry at a fixed temperature. Li and Wu 16J 
modeled the ordinary and extraorinary refractive indices 
wavelength dependence based on the extended Cauchy for- 
mulae. Their model is described by the following equa- 
tions: 

Bp Cf , . 

n e =A e -^ + ^, (23) 

_ A B o C 

n Q - a - + -jj. 

The coefficients appearing in equations (|23p and 
were obtained [16] by fitting experimental data. For 5CB 
at 25.1 °C they are given in the tables below. 



(24) 







C e 


1.6795 


0.0048 [im 1 


0.0027 ^m 4 




A 


B 


Co 


1.5174 


0.0022 [im 1 


0.0011 ^m 4 



In figure 5 we show the ratio a = n /n e as a function 
of the wavelength as obtained from equations (|23|) and 



Since both temperature and wavelength cause a to 
change we can summarize their effect on the light paths by 
studying the geodesies for different values of a. Substitut- 
ing the metric @ in © and this one in ([5]) we can calcu- 
late the geodesies for different values of a. As described in 
Section II, the geodesic equations ((8|) have exact solutions 
for the k = 1 case. The remaining cases can be solved by 
a numerical method. In figures 6 and 7 we show the ef- 
fects of the variation of the parameter a on the light paths 
near the k = 1 defects, using the exact solution of Section 
II. In figure 8 we show the same effects for the k = — 1 



disclination, using the Runge-Kutta numerical method to 
solve the geodesic equation. In all cases, the solid line cor- 
responds to a=0,8912, the dotted line to a=0,9120 and 
the dash-dotted line to a=0,9355. These values, for 5CB, 
probed by a 589 nm light beam, correspond to the tem- 
peratures of 290 K, 300 K and 305 K, respectively. Notice 
that as a approaches 1 the light paths straighten out, as 
it should. 




Fig. 6. Influence of the parameter a on the light trajectories 
in a nematic liquid crystal with a disclination k = 1 and c = 0. 
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Fig. 7. Influence of the parameter a on the light trajectories in 
a nematic liquid crystal with a disclination k = 1 and c = n/2. 
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Fig. 8. Influence of the parameter a in the light trajectories 
in a nematic liquid crystal with a disclination k = — 1 and 
c = 7T/2. 



4 Conclusion 



Topological defects in nematics cause light passing by to 
deflect, as shown in [6|. In that article, we associated the 
light paths to geodesies in a curved space specified by the 
defect. The deflection is due to the particular orientation 
of the director field associated to the defect, which may be 
translated into curvature. The intensity of the deflection 
depends on the ratio a between the ordinary and extraor- 
dinary refractive indices, which, in turn, depend on the 
temperature of the liquid crystal and on the wavelength 
of the light. Taking as example 5CB, which has been ex- 
tensively characterized with respect to temperature and 
wavelength dependence of the refractive indices p^4lfl6] , we 
solved the geodesic equations for a realistic range of val- 
ues of a corresponding to temperature and/or wavelength 
variation. The graphical result illustrates the influence of 
these parameters on the light deflection caused by the de- 
fect. The further a gets from 1 the strong is the deflection. 
This can be achieved by either lowering the temperature 
or shortening the wavelength. In conclusion, the study of 
the influence of measurable physical parameters, like tem- 
perature and wavelength, helps us to understand better 



the behavior of light propagation in liquid crystals where 
topological defects are relevant. 
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